Short-time quantum detection: probing quantum fluctuations 
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In this work we study the information provided by a detector click on the state of an initially 
excited two level system. By computing the time evolution of the corresponding conditioned prob- 
ability beyond the rotating wave approximation, as needed for short time analysis, we show that 
a click in the detector is related with the decay of the source only for long times of interaction. 
For short times, non-rotating wave approximation effects, like self-excitations of the detector, for- 
bid a naive interpretation of the detector readings. These effects might appear in circuit QED 
experiments. 
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Quantum detection theory was created to study the 
behavior of detectors in presence of radiation [1 . Highly 
satisfactory up to date, it relies on the conspicuous ro- 
tating wave approximation (RWA), which neglects the 
so-called counterrotating terms, irrelevant in most cases. 
These terms give important contributions for strong 
atom-field couplings and very short times as compared 
to the system time scale, meaning that for any effect be- 
yond RWA to be acknowledged our measurements must 
be very precise and the observation times quite fast. This 
is particularly problematic for Quantum Optics experi- 
ments, due to the very small matter-radiation coupling 
and the fact that observation times must be at the fem- 
tosecond scale for most cases (nanosecond for hyperfine 
qubits), which is ridiculously small for current experi- 
ments (~ jis for trapped ions [2]). 

However, circuit QED [3] provides a framework in 
which those phenomena are accessible to explore. By 
using superconducting qubits as artificial atoms coupled 
to a transmission line, one gets a system which be- 
haves analogously as a 1-D radiation-matter interaction 
in Quantum Optics, while working at the microwave fre- 
quency range [4 . In this setup, parameters can be easily 
tuned, and the coupling between qubits and transmission 
line can be modulated up to ultrastrong levels. [5] ^6 Be- 
sides, fast qubit state readout (~ ns) is possible using a 
pulsed DC-SQUID scheme ^ . 

In those conditions, phenomena beyond RWA have al- 
ready been reported [8], [9]. In this regime Glauber's 
theory is no longer valid and quantum detectors should 
be described by a non-RWA model. 

A counterintuitive direct consequence of the break- 
down of the RWA is that a detector in its ground state 
interacting with the vacuum of the field has a certain 
probability of getting excited and emitting a photon. 
There is however not a real consensus on the physical 
reality of this effect. Introducing counterrotating terms 
is interpreted by some to be a problem as the processes 
described by those terms seem virtual. It seems coun- 
terintuitive to accept that a detector in a ground state 
in the vacuum could get excited. As a matter of fact. 



there have been attempts of suggesting effective detector 
models by imposing this phenomenon to be impossible 

We should however point out here that there is no 
problem whatsoever with energy conservation, as unitary 
evolution of the states under a hermitian hamiltonian di- 
rectly guarantees it. Any peculiar effect at that respect 
is linked to the fact that the initial state considered has 
not definite energy, since the state "detector and field 
in their ground states" is not an eigenstate of the full 
Hamiltonian beyond RWA. 

In this work we would rather step out of this discussion 
and by taking the theoretical models coming from circuit 
QED, without imposing any additional constraints, we 
will study the following setup: a source S initially ex- 
cited, a detector D initially in the ground state and both 
interacting with the electromagnetic field in its vacuum 
state. If the detector clicks at a given time, does it mean 
that the source is now in the ground state? This prob- 
lem amounts to compute the probability of decay of the 
source, conditioned to the excitation of the detector. We 
will show that, unlike Glauber's RWA detector in which 
this conditioned probability would be equal to 1 at any 
time, this circuit QED detector only achieves this value 
at long times, due to the impact of non-RWA effects. 

More precisely, our model has two superconducting 
qubits, S and with two levels g and e and separated a 
distance r. Let us consider an initial moment t = where 
S is excited, D (which represents the detector) is in its 
ground state and there are no excitations in the trans- 
mission line, which will be open, so enabling a contin- 
uum of modes. Representing the states in terms of qubit 
(6*, D) and field (F) free eigenstates with the notation 
= \SDF)^ the initial state would be |^)^^o ~ \^9^)' 

After a certain time t, if we measure qubit D and it 
results excited, that would naively lead us to think S 
has decayed and produced a photon which has then later 
been absorbed by D. We intend to proof otherwise by 
quantifying what information about the state of S can be 
extracted by knowing qubit D state after a certain time 
t. For that we will compute the probability Vsg/oS'^) 
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S to have decayed at a certain instant t, conditioned we 
have measured D excited at that same moment: 

V rf^-^bH _ EFl(geF|e-'^*/^|egO)|^ 

^[^e*] being the probabihty of 5* being in the ground state 
and D excited and 7^[*e*] the total probabihty of excita- 
tion of D. 

From here on, we will consider the following Hamilto- 
nian [ill [HI: 
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A={S,D} 

Hi = - dAV{xA)cjt 

A^{S,D} 



(2) 



Here xa corresponds to the position of the qubit A, HVLa 
is the gap between levels for qubit A and V refers to the 
1-dimensional field which expands as: 



V{x) = i / 

J — ( 



dk^jNujke^'^'^ak 



H.c, 



(3) 



This field has a continuum of Fock operators [a/c,a^/] = 
S{k — k')^ and a linear spectrum, uok = v\k\^ where v is 
the propagation velocity of the field. The normalization 
and the speed of photons, v = (c/)~^/^, depend on the 
microscopic details such as the capacitance and induc- 
tance per unit length, c and I. Note that this model, a 
dimensional reduction of the matter- radiation hamilto- 
nian with two-level atoms and the electromagnetic field, 
is formally equivalent to the Unruh- de Witt model. [13]. 
For our calculations, we will make use of the interaction 
picture, so we let the initial state |e^0) evolve for a lapse 
of time t as: 

m)) = Ui{t) \egO) = T{e-^f^dt'H,it^)/n^ 

= / \egO) + X \geO) + ^ Ai,^ \gglk) + Yl ^^'^ 1^^^^) 

k k 

+ \eg2kk') + \ge2kk') + (4) 
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Here and in the following we will only make explicit the 
terms that contain contributions for the probabilities up 
to d\. For example, terms with 3 or more photons in the 
amplitude will be excluded, as they give contributions of 

Having A4{t;nF) = {neF\tlj{t))^ the probabilities 
needed for the computation of Vsg/D^it) ^ can be writ- 



ten down using Q as: 

^be*] = E \{9eFmt)\egO)\^ = ^ \M{t;gF)\^ 

F F 

= \Xf+J2\B2\' + ... (5) 
^[*e*] = J2 \{neF\Ui{t)\egO)f = ^ \M{t;nF)f 

n,F n,F 

= |X|2 + ^|Bip + ^|i?2|' + ... (6) 
The first building block needed is Note that: 



V^geo] = \{geO\Ui{t)\egO)\^ = \M{t;gO)(^ 



\X\ 



(7) 



To evaluate |Xp up to fourth order in perturbation the- 
ory, one must consider X has no contributions neither 
for orders or 1, so the calculation must be performed 
for orders 2 and above. As a fact, order 2 alone will 
suffice. This calculation has been already performed in 
the appendix of [14 , where the dimensionless coupling 



Ka 



Ad\N 



with A = {S^D}, acts as the perturbative 



parameter. We will restrict to times where i^^cj^t <C 1, 
where our perturbative approach remains valid . 

In Fig. [l] we sketch the evolution of the probability 
V[geo] with time, and its dependence with the coupling 
and the distance between qubits. Typical values for cou- 
plings and distances for a setup in circuit QED are con- 
sidered from here on. At these early stages V[geo] is highly 
oscillatory in time. For a given time, the probability al- 
ways grows with the coupling strength but depends of the 
distance in different ways. To proceed with the calcula- 
tion of V[ge^] , the term B2 must be evaluated. As the final 
bare state associated with this term has two photons, this 
implies automatically that orders and 1 are discarded. 
Once again, order 2 alone fits. Once calculated, it must 
be squared and summed, splitting into two terms, a "di- 
rect" one, just the product of the square of the emission 
amplitudes, and a "crossed" one, which looks like a pho- 
ton exchange. The summation of the direct terms im- 
plies the appearance of expected divergences which can 
be resolved using a regularization procedure as the one 
sketched in the appendix of [15] . This procedure requires 
the times of analysis to be larger than a certain cutoff 
time to, which in this case is related with the typical size 
of a superconducting qubit d 10~^ m [16] and the prop- 
agation velocity of the field quanta: v ^ 10^ m/s. Thus, 
to = d/c 1 = 10~^^ s far below the times considered 
in this work. A detailed treatment of the procedure re- 
lated with the emission probabilities at short times will 
be published elsewhere. 

Notice that B2 is only non-zero beyond the RWA. 
In Fig. I2] we compare V[ge^] with V[geo] and the im- 
pact of tms non-RWA contribution is seen in the sub- 
nanoseconds regime for a large value of the coupling 
strength. At larger times, the impact diminishes and 
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FIG. 1. a) P[egO] in front of t in s for three different values 
of the distance between qubits 27Vj = 0.1 (dotted, squares, 
blue), 0.3 (dashed, crosses, green) and 0.5 (solid, circles, 
black). For all cases the coupling strength is Ks = Kd = 
7.5 • 10~^ and the qubit frequency Q/(27r) = 1 Ghz. b) P[ego] 
in front of t in s for three different values of the coupling 
strength K = Ks = Kd = 6 - 10"^ (dotted, squares, blue), 
7.5 • 10"^ (dashed, crosses, green) and 9 • 10 ^ (solid, circles, 
black) . For all cases 27r^ = 0.3 and Q/{27v) = 1 GHz. 
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FIG. 2. a) P[eg*] (solid, blue, circles) and P[ego] (dashed, 
green, crosses) in front of t in s with a distance 27r J = 0.5, a 
coupling strength K = Ks = Kd = 3- 10~^ and a qubit fre- 
quency of Q/(27r) = IGHz {Q = Qs = ^d). The difference 
between the two graphs is the non-RWA term 5^ | 1 ^ • 



The last probability of interest, 7^[*e*]7 needs XI l^ip, 
which is again a completely non-RWA contribution. For 
that case the situation gets more complicated, as there 
are interfering processes of orders 1 and 3 leading to that 
final state. This calculation has been the focus of a recent 
work [17] . The four diagrams contributing to ^ 1 5i p up 
to fourth order in perturbation theory can be seen in fig 



(|3|. The leading order contribution is just the probability 
of self-excitation of the detector, (fourth diagram in fig. 
|3| and the other contributions come from the interference 
of this diagram with the other three. In particular, in- 
terference with the third diagram is crucial for causality 
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FIG. 3. The diagrams contributing to X, Bi and B2. X 
represents the amplitude for photon exchange between source 
and detector, while B2 is just the amplitude for two single 
photon emissions, one at each qubit. The leading order con- 
tribution to Bi is the amplitude for a single photon emission 
at the detector qubit (fourth diagram), but third-order one- 
loop corrections (first two) and a photon exchange accompa- 
nied by an emission at the source have to be also taken into 
account. Bi and B2 are completely non-RWA diagrams. 

With the previous probabilities computed we can 
finally address the conditioned probability Vsg/D^i'^)^ 
which can be calculated as shown in ([T]). Note that 
in the RWA Vsg/oS^) = 1 at any time, since 7^[*e*] = 
^[^e*] = ^[geo]' We have seen however that non-RWA 
contributions to 7^[>^e*] have a sizable impact in the sub- 
nanosecond regime. The effect of these contributions to 
the evolution of Vsg/oS^) ^een in Figs. [4]and[5j 

where the consequences of changing the coupling and tne 
distance between qubits are considered. The first thing 
we notice in fig. [4] is that for short times the informa- 
tion provided by the detector is not very much related 
to the state of the source, that is, self-excitations and 
other non-RWA phenomena dominate over the photon 
exchange between source and detector. For the cases 
considered, only at interaction times t ^ Ins the 
conditioned probability converges to the RWA prediction, 
that is, the excitation of the detector is a reliable way to 
detect the decay of the source. Since the non-RWA con- 
tributions are more relevant for large couplings and short 
distances, the convergence is faster as the distance grows 
and the couplings diminish, as can be seen in Figs. |4]and 

m 

The above theoretical results could have an impact in 
real experiments of circuit QED. In particular, a typical 
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FIG. 4. Vsg/Dei^) front of Kt for three different values 
of the coupUng strength oi K = Ks — Kd — 7.5 • 10~^ 
(sohd, blue, circles), 1.5 10~^ (dashed, green, crosses), 7.5- 
10~^ (dashed, black, squares). In the three cases 2ti^ — 1 
and VL/{2t^) = 1 GHz (Vl^VLs^ ^d). 
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FIG. 5. Vsg/Deif) front of t in s for three different 
values of the distance 27r^ = 0.5 (dotted, blue, circles), 0.75 
(dashed, green, crosses), 1 (solid, black, squares). In the three 
cases the coupling strength K — Ks — Kd = 1.5 • 10~^ and 
Q/(27r) 1 GHz = Qs = ^d). 



setup to measure the internal state of a flux qubit cou- 
pled to a transmission line consists of a SQUID surround- 
ing the qubit. Although the total measurement process 
could take up to tens of nanoseconds, most of the time 
the coupling SQUID-qubit is much stronger than K, [7] 
and the dynamics qubit-transmission line is effectively 
frozen. Thus this dynamics is only important during the 
activation of the SQUID, a process that may be in the 
nanosecond regime. Therefore, if at t = the system is 
in the initial state of \egO) (see [17j) for the preparation 
of such a state) and then we measure the state of the 



qubit, our results show that a click of the detector could 
be as related to a self-excitation as to an absorption of a 
photon emitted by the source. 

To conclude, we have considered a system of two su- 
perconducting qubits coupled to an open transmission 
line, which can be suitably described in the framework of 
1-D QED with two-level (artificial) atoms. By using per- 
turbation theory in the interaction picture we have com- 
puted at short times the probability of excitation of a de- 
tector qubit in the presence of an originally excited source 
qubit subjected to a possible decay. Our main result is 
that for typical circuit QED parameters, a significative 
amount of time is needed before we start trusting the de- 
tector results as informative with respect to the source. 
This is an effect of the breakdown of the RWA in circuit 
QED, which gives rise to non-RWA phenomena like self- 
excitations. By neglecting the counterrotating terms, the 
sophisticated reality of the short-time phenomena stays 
hidden and, in particular, a total reliability on the in- 
formation coming out of the detector is wrongly derived 
for all time-scales. It is noteworthy to mention that this 
result applies in general to other setups and detectors, 
although for extremely short times, and it is in the case 
of circuit QED where it might affect the interpretation 
of coming experimental results. 
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